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Abstract 

The theoretical description and the properties of Frenkel excitons in non-covalently bonded molecular 
aggregates are reviewed from a multi-exciton perspective of dissipative quantum dynamics. First, 
the photophysical and quantum chemical characterization of the monomeric dye building blocks is 
discussed, including the important aspect of electron-vibrational coupling within the Huang-Rhys 
model. Supplementing the model by the Coulombic interactions between monomers, the description 
of aggregates in terms of excitonic or vibrational-excitonic bands follows. Besides of giving rise to 
complex absorption and emission line shapes, exciton-vibrational interaction is responsible for en- 
ergy and phase relaxation and thereby limits the size of coherent excitations in larger aggregates. 
Throughout, emphasis is put on the electronic three-level model as a minimum requirement to de- 
scribe nonlinear spectroscopies including effects of two-exciton states such as excited state absorption 
and exciton-exciton annihilation. The experimentally observed characteristics of stationary absorp- 
tion and fluorescence spectra of aggregates as well as their temperature dependence are discussed. 
Examples for ultrafast spectroscopic experiments including pump-probe studies, photon echo and two- 
dimensional spectroscopy are presented and results on the size of coherence domains and on intra- 
and interband relaxation are given. Finally, experimental signatures for exciton-exciton annihilation 
and their analysis with respect to the mobility of excitons are described. 

Key words: photophysics, Frenkel excitons, nonlinear spectroscopy, electron-vibrational coupling, 
open system dynamics 



1. Introduction 



Organic solar cells are a highly active research field since they hold the promise to provide a cheap 
source for renewable energy. In such cells the conversion process of a photon into electric energy 



consists of several steps (Brabec et al. (20011; Gregg (2003|; [Ayzner et al. (2009)). First the photon 
is absorbed by dye molecules or a polymer and an exciton is generated. The exciton has to reach 
an interface between two materials with shifted band gaps. Here charge separation can take place. 
The generated charges have then to migrate to the contacts. In commonly used organic materials 



the exciton diffusion length is limited to about 10 nm (Scully and McGehee (2006); Mikhnenko et 
al. (2008)). This short distance poses severe restrictions to the design of organic solar cells. A wide 
spread concept to deal with the short exciton diffusion length arc bulk hctcrojunctions in which two 



materials form strongly interleaving domains (Brabec et al. (2001)). However, not too many materials 
are suitable for bulk heterojunctions. Another option is to design organic materials with a high 
degree of order which provide a long exciton diffusion length. For this purpose aggregates are highly 
attractive since, in addition to high exciton mobilities, they exhibit also strong absorption bands in 
the visible spectral region (Mobius (19951; Wurthner et al. (2011)). Aggregates are highly ordered 
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Figure 1: (a) Different processes in aggregates composed of three level monomers, (b) Excitation and internal conversion 
process in a model of displaced harmonic oscillators describing the nuclear motions of the monomer. Notice that the 
internal conversion process usually proceeds via intermediate electronic states as indicated by the shaded area. 



supramolecular structures of dye molecules which form by self-assembling. The monomers are bound 
by van-der-Waals forces and often in addition by hydrogen bonds. To understand and design their 
cxcitonic properties is crucial for their application in organic solar cells as well as in other photonic 
devices. The theoretical description and spectroscopic investigation of these excitonic properties are 
the topic of this review. 

Molecular aggregates store and transfer electronic exciton energy by means of Frenkel excitons 
(Knoester and Agranovich (2003); Beljonne et al. (2009); May and Kiihn (2011)). Frenkel excitons 
are electron-hole pairs localized at the same molecular unit (monomer) of an aggregate. They have 
to be distinguished from Wannier-Mott excitons common in solids, where the electron-hole separation 
reaches hundreds of nanometers. Further, so-called charge transfer excitons play a role in molecular 
systems, provided that a means for charge separation e.g. between parts of a larger molecular unit 
is given. Exciton formation has been observed first by Jelly and Scheibe in 1936, who detected a 
narrowing of the relatively broad absorption band of pseudo-isocyanine (PIC) upon increasing its 
concentration in solution (for reviews, see Mobius (1995) and Wurthner et al. (2011)). Nowadays 



it is recognized that the concept of excitons is not restricted to molecular aggregates, but is a basic 
mechanism for electronic excitation energy transfer (EET) in various nanoscale and biological systems 
( |Scholes and Rumbles| ( |2006[ )). 

Traditionally the theoretical description of electronic EET has been based on the model of elec- 
tronic two- level systems, coupled via dipole-dipole interactions (Forster (1948); Kasha et al. (1965)). 
Considering the aggregate to be build up of monomers which retain their chemical identity, the dipole- 
dipole coupling is an approximation to the full Coulomb interaction between the electron densities of 
the monomers. Its limitation is obvious and especially severe in case of tightly packed aggregates. In 
recent years the advances in computational power and numerical methods have made it possible to 
determine the full Coulomb coupling between monomers from first principles, even for rather large 
dye molecules. 

Scrutinizing the validity of the two- level approximation is equally important. In fact for the 
understanding of linear spectroscopies this model is well suited. However, for the description of 
nonlinear spectroscopy as well as of processes triggered by higher laser intensities the incorporation 
of higher excited electronic states might become vital. Here, more photons are absorbed, usually at 
similar energies. Viewed from the monomeric perspective of Fig. [TJd, after, e.g. a So — > Si excitation, 
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a second photon may excite a Si — > S„ transition, where n could be of the order of 10 or higher. 
Thus, strictly speaking the modeling of all processes following this excited state absorption (ESA) 
would require the use of n electronic levels for each monomer. However, according to Kasha's rule it 
can be expected that the S n state population quickly relaxes down to the Si state. Thus to a good 
approximation the monomer can be described by an electronic three-level model, supplemented by a 
rapid internal conversion mechanism between the S„ and Si states. 

There is a second type of excitation, which can be reached upon absorption of two photons, i.e. 
two monomers might be excited to the Si state yielding a nonlocal doubly excited (two-exciton) 
state of the aggregate. Besides the importance of these local and nonlocal doubly excited states for 
spectroscopy, their interplay has some consequences for the dynamics of EET as illustrated in Fig. [l^. 
Consider a nonlocal double excitation (upper panel) . The two excitations can transfer independently 
through the system, but owing to the Pauli principle a monomer Si state cannot be doubly excited. 
However, the two excitations can combine to yield a local S„ excitation (exciton fusion, middle panel). 
The life time of the S„ state, however, is usually short and internal conversion will quickly lead the 
system back to the Si state. Hence, in total one electronic excitation has been dissipated into heat, i.e. 
vibrational motions. This so-called exciton-exciton annihilation (EEA) mechanism is an important 
limiting factor for device applications as it restricts the amount of excitation which can be taken up 
and transported by the aggregate. 

The microscopic origin of EEA is the breakdown of the Born-Oppenheimer approximation, leading 
to non-adiabatic couplings between the adiabatic electronic states triggered by nuclear motions. How- 
ever, even within the Born-Oppenheimer approximation there will be effects of electron-vibrational 
coupling on the dynamics and spectroscopy of Frenkel excitons. This includes the appearance of vi- 
brational progressions in particular for isolated higher frequency vibrational modes in the absorption 
and emission spectra, which combined with the Coulomb coupling may result in rather complex line 
shapes. The coupling to the dense manifold of low-frequency vibrational modes will cause an overall 
line broadening, where one can distinguish between static and dynamic disorder according to the time 
scale of nuclear motions with respect to relevant times of the exciton motion. In the time domain the 
effect of exciton-vibrational coupling (EVC) is to cause exciton coherence dephasing and population 
relaxation. Thus EVC limits the size over which a coherent excitation can be established within an 
aggregate and therefore is crucial for exciton transport. 

In this review we will discuss exciton dynamics and spectroscopy focussing on some of the points 
just raised, namely the three-level nature of the monomeric building blocks, exciton-exciton annihi- 
lation, and exciton-vibrational coupling. We will restrict the discussion to molecular aggregates; for 
an account on excitons in biological pigment-protein complexes and artificial nanoscale systems like 
polymers we refer the reader to other contributions in this book. Furthermore, we will not treat 
the issue of exciton dissociation and free charge carrier formation which is crucial for the function of 



excitonic organic solar cells (Gregg (2003)); for a review on recent theoretical advances see Burghardt 
et al. ( 2009 ) . After a brief account on various aspects of Frenkel exciton theory in Section 2j spectro- 



scopic applications will be discussed in Section [3j which highlight the unique properties of excitons in 
molecular aggregates. 



2. Theory 



2.1. Molecular Aggregate Hamiltonian 

Separating the aggregate into its constituent monomers, which are coupled by the Coulomb inter- 
action, the aggregate Hamiltonian becomes (Renger et al. (2001); [May and Kiihn (2011)) 
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Here, \(p ma ) denotes the ath adiabatic electronic state of the mth monomer, which depends on the 
nuclear coordinates R m = (R m .i,R m ,2) . . .) of the monomer, i.e. ip ma — (p ma {Rm)- This dependence 
is found also for the matrix elements entering Eq. ([I]). For the monomer contribution we have 



H„ 



-^m,a(-^m) -^m 4~ ^)nfl(^m) 



(2) 



Here, T m and V ma (R m ) denote the nuclear kinetic energy and potential energy surface, respectively. 
Further, in Eq. ([T]) m . a b is the operator of non-adiabatic coupling between different adiabatic 
electronic states and it is assumed that & m .aa — . The second sum in Eq. (JlJ contains the Coulomb 
coupling between monomers m and n. It can be expressed in terms of transition densities 
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as follows 



J mn (ab, cd) 
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In Eq. (|3| g^\x) is the electronic transition density for a transition between the adiabatic states 
\tpma) and \(p m b) at monomer m and the are the atomic numbers. Often the Coulomb coupling is 
treated within the dipole approximation assuming that the spatial extension of the transition density 



is small as compared with the distance between the monomers ( |May and Kiihn (2011 1). This gives 



J^(ab,cd) = 



be 



IX I 3 



(5) 



Here, X m „ is the distance vector between the considered monomers (e.g. center of mass) and the 
d m at, are transition dipole moments. 

Considering the influence of nuclear motions on the EET dynamics one has to distinguish between 
the effect on the intramolecular transition energies (via V ma {R m )) and on the Coulomb coupling. 
For aggregates in solution there will be an additional dependence of both quantities on the solvent 
coordinates. As discussed in May and Kiihn (2011) there are different models to cope with this 
situation (cf. Sec. 2.4). This includes the description in terms of combined local intramolecular 



vibrations and intermolecular as well as solvent normal modes, the use of global normal modes for the 
total aggregate plus solvent system, and the classical description of nuclear motions. While the former 
two cases have been in use for quite some time, explicitly accounting for the thermal fluctuations by 



virtue of molecular dynamics simulation of aggregates has become feasible only recently (see, e.g. Zhu 
eF^([2007||2008bl). 



2. 2. Frenkel Exciton Hamiltonian for Coupled Three-Level Monomers 

EET is commonly discussed on the basis of an electronic two-level approximation for the monomers. 
While this is sufficient for the discussion of linear optical properties, the understanding of nonlinear 
optical experiments as well as of the dynamics at higher excitation densities requires to take into ac- 
count the fact that organic dye molecules forming the aggregate show ESA at wavelengths comparable 
to that of the Sq — ¥ Si transition. A minimal model which can mimic this situation is an electronic 
three-level description (cf. F ig. [Tjp). The first stud y of optical properties of linear three-level aggre- 
gates has been presented by Knoester and Spano (19951. These authors showed that the coupling 



between local and nonlocal double excitations results in distinct features of the two-photon absorption 
spectrum. A region of coupling parameters has been identified which will lead to a strong distortion of 
the pump-probe spectrum. In fact under certain conditions (harmonic limit) the third-order response 
will vanish all together (Knoester and Spano (1996)). This issue has subsequently been addressed for 
arbitrary aggregate geometries on the basis of the anharmonic exciton oscillator picture (Kiihn et al. 
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( 1996b l). In particular a time-domain four- wave mixing technique has been proposed to study the 
interplay between local and nonlocal double excitations (Kiihn and Mukamel (1997)). 

In the following we will specify the aggregate Hamiltonian, Eq. (|lj), to the case of coupled three- 
level monomers. Restricting the summation indices to a = (g, e, /) leads to a site-diagonal term which 
contains the non-adiabaticity operator mixing all states. However, in practice the non-adiabatic 
coupling to the ground state can be neglected, i.e. Q m ,ga ~ (a — (e, /)). Concerning the coupling 
Q m je the restrictions of the three-level description becomes apparent (cf. Fig. [TJ}). Normally, the 
state \<Pmf) will not be the second excited state, but a certain higher excited state, S„, which happens 
to have about the same transition frequency with respect to Si as the So — > Si transition. In other 
words, there will be a sizable number of electronic states, which are energetically in between \(p me ) 
and \tpmf) and it is rather likely that the electronic relaxation proceeds via these states. Since not 
much is known about these states and their couplings, Q m j e has to be viewed as an effective coupling. 

Next we discuss the Coulomb integrals J mn (ab,cd), which can be classified as follows (May and 
Kiihn (2011)): First, there are interactions between charges (J mn (aa,aa), J mn (ab,ab)) and between 
charges and transitions (J mn (ab,ca) or J mn (ab,bd)). They contain the effect of aggregation on the 
adiabatic energy levels. If no first principles determination of the aggregate Hamiltonian is attempted, 
this effect is usually considered as being included in the definition of the adiabatic states. Note that 
these terms are also held responsible for the exciton-exciton interaction which can lead to the formation 



of bi-excitons (Spano et al. (1991); Juzeliunas and Reineker (1998); Knoester and Agranovich (2003)) 



Second, there are matrix elements of the type J mn (aa, bb) that present interactions between simultane- 
ous excitations and de-excitations of both monomers. From the energetic point of view this will lead to 
a rather off-resonant process. Neglecting these terms is usually called Heitler-London approximation. 
Within this approximation the number of excitations will be conserved. Finally, matrix elements of 
the type J mn (ab,ab) and J mn {ab,cc) which give the interaction between excitation and de-excitation 
of the monomers. In contrast to the previous case this interaction facilitates resonant EET and there- 
fore it is often the only term that is considered in modeling of EET. For the electronic three-level 



J mn (eg,eg) = ./, 



(eg) 



system there are two types of matrix elements which need to be considered, i.e. 
which facilitates resonant EET involving So — >Si transitions and J mn (fg, ee) = Jmn which involves 
simultaneous So — >Si and Si — >S„ transitions. 

Using this approximation the aggregate Hamiltonian, Eq. ([!]), can be recast into a simpler form. 
To this end a classification of the electronic states according to the number of excitations which are 
present, e.g., as a consequence of the interaction with the light field is introduced. For the description 
of third-order nonlinear optical spectroscopy it suffices to include the ground state 



n 



<Pr. 



(6) 



the state carrying a single excitation 

|m) = \<p me ) Yl Wng) , (7) 
and the one carrying two excitations which can be either local or nonlocal 

\mn) = (1 - 5 mn )\ip me )\ip ne ) Yl \fkg) + <5 > n \^ g ) ■ (8) 

k^m,n ky^m 

Introduction of this excitation expansion allows to rewrite the Hamiltonian as follows 

H = H<®+H<»+H&>+Hl 1 -*K (9) 
Here the ground state Hamiltonian contains just the contributions from nuclear motions 

H<® = ^H m , g \0)(0\=£o\0)(0\. (10) 
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For the single excitation manifold we have 



H {1) = 5> m „(£o + U m , eg ) + &\m)(n\ , (11) 

where 

U m ,ab — V m , a — Vm,h (12) 

is the nuclear coordinate dependent electronic gap coordinate, i.e. U m a b(R m ) Finally, we have for the 
double excitation manifold 

fji 2 ) — ^ ^ {<5 m fc<5„/[(l — S mn )(£() + U„ lyeg + Un^eg) + 8 mn (£o + U m Jg)] 
l>k n>m 

+ (1 - <W)(1 - 5 kl )[J ( Z ] 5m + 4* 9) 5 mk + 4 e n 9) S lm + Jl e n 9) 5 kn ] 

+ S mn (l - S u )[4 f ^5i m + J^hk^kl)^ . 
The non-adiabatic coupling has only on-site terms, i.e. 

H (i-z) = J2&mje\mm){m\+h.c. (13) 

n i 

The same holds true for the electric dipole coupling to the laser field which reads 

#fieid(*) = -E(i)^d m , eff |m)(0| -E(t) ^{5 mn d mJe \mm)(m\ 

m n>m 

+ (1 - 5 mn )[d m . eg \mn)(n\ + d n , eg \mri)(m\]} + h.c. . (14) 

Eqs. ( 10p4 ) provide a complete description of the laser-driven exciton dynamics up to doubly excited 
states. By virtue of the electronic gap coordinates, Eq. (12), and the non-adiabatic coupling the 
Hamiltonian still depends on the nuclear degrees of freedom. This will lead to energy and phase 
relaxation processes within the excitation subspaces (intraband) as well as to transitions between them 
(interband). Leaving aside the non-adiabatic transitions and considering cases where the Coulomb 
coupling dominates the effect of the EVC, it is customary to introduce Frenkel exciton states by 
separate diagonalization of the single (JV = 1) and double (N — 2) excitation Hamiltonians for a fixed 
nuclear configuration (i.e. at the ground state equilibrium) 

HW\a N ) =E aN \a N ). (15) 

The exciton states can be decomposed into the basis of the local excitation states as follows: 



a. 



i)=Y^ C ^\m), (16) 



|«2> = ^2 C mn ^ 2 \mn) . (17) 

Notice that if -/V mo i is the number of monomers, one has N mo \ one-exciton and iV mo i(iV mo i + l)/2 
two-exciton states which are often called bands. 

Using this expansion, the non-adiabatic coupling, Eq. (13), and the aggregate-field interaction, 
Eq. (14), can expressed as couplings between different exciton manifolds. 

In general, Eq. (15), has to be solved numerically. However, for special geometries analytical 



solutions are available. This includes the dimer, the linear open and closed chain (May and Kiihn 



G 



( 201 1 ) ) and certain tubular aggregates ( Knoester ( 2006 ) ) . For example, for a homodimer with identical 
site energies Hu} eg and coupling J, the cigencnergies of the one-exciton manifold are given by 

E ai=± = ^±J (18) 

and the eigenstates are 

|a 1 =±) = i(|l)±e-^W|2». (19) 

For the case of a linear aggregate of N mo ] monomers having identical So —*Si transition energies huj eg 
and nearest-neighbor coupling J one has 

E ai = hw eg + 2Jcos( AT na \ - ) ai = l,...,N mo i (20) 



N mo i + 1 



with the eigenstate expansion coefficients 



Cm -=v^^ s H m £^) (21) 

Notice that linear aggregates are often considered by using periodic boundary conditions, which enables 
one to introduce the eigenstates as Bloch states having quasi- momentum k (Knoester ( 2006| >). 



2.3. Quantum Chemical Methods 

There are several strategies for the determination of the parameters of the Frenkel exciton Hamilto- 
nian. On an empirical level, fitting to experimental data can be attempted. With increasing computer 
resources and the development of sophisticated methods for the solution of the electronic Schrodingcr 
equation even for larger chromophorcs in solution, a first principles based Frenkel exciton Hamiltonian 



can be constructed. Although semiempirical methods still enjoy popularity (Fron et al. (2008)), in 



particular time-dependent density functional theory (TDDFT) plays a prominent role as it often gives 



a reasonable compromise between accuracy and numerical efficiency (for a critical overview see Liu 



et al. (2011)). Here one can distinguish between two approaches: First, the aggregate is treated as 
a supermolecule, usually restricted to a representative dimer. Based on the splitting between excited 
states and their dependence on the dimer geometry, coupling parameters can be deduced. 

Second, the aggregate is viewed as being composed of monomers whose monomeric electron den- 
sities are interacting via the Coulomb potential. Thus one follows the line of the Frenkel Hamiltonian 
to combine monomer transition energies and dipoles with calculations of Coulomb matrix elements. 
There are several methods for obtaining J mn (ab, cal) in Eq. Q. Scholes and coworkers (Krueger et al. 



(1998)) have been the first to propose a practical method for a first principles based calculation of 
the Coulomb coupling. Their TDC (transition density cube) approach coarse-grains the space for in- 
tegration of the transition densities, thus expressing the interaction as a sum over pairs of TDCs. An 



alternative has been suggested by Renger and coworkers ( jMadjet et al. (20061). Their TrEsp (tran 



sition charge from electrostatic potential) method employs atomic partial charges in the summation 
of the Coulomb interaction. These partial charges are obtained by fitting the electrostatic potential 
associated with the monomeric transition densities. If q^"" (a, d) denotes the Ith charge associated 
with the transition density the Coulomb integral becomes 

|K/-Kj| 

where R/ is the position of the partial charge. As compared to the TDC method this approach has 
the advantage that essentially converged results can be obtained using some tens of partial charges. 
More recently direct implementations of the calculation of the Coulomb coupling became available, 
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which utilize the integration and pre-screening tools available in quantum chemistry packages ( Fiickel 
eTaI1p008| ). 

Studying the higher excited states which are relevant for two photon transitions, one faces a serious 
computational challenge. The Si — > S„ transitions often involve transitions between molecular orbitals 



with double excitation character (Pabst and Kohn (2008)) and therefore cannot be treated by TDDFT 
methods. Recently, we have shown that the MRCI/DFT (multi-reference configuration interaction) 
approach developed by Grimme and Waletzke ( 19991) can provide guidance for the assignment of Si — > 
S„ transitions in large chromophores such as perylene bisimides (Ambrosek et al. (2011)). 



2.4- Exciton- Vibrational Coupling 

The coupling between electronic and nuclear degrees of freedom leads to the appearance of the 
energy gap coordinate U mta b, Eq. (12), in the monomer Hamiltonian. Furthermore, upon aggregation 



the Coulomb coupling between the monomers depends on the intermonomer distance and its variation, 
e.g., due to thermal fluctuations, will have an influence on the exciton dynamics. Besides theses intra- 
aggregate coordinates, there will be an influence on the degrees of freedom of the solvent environment. 
There are various strategies for incorporating the effect of EVC. Of considerable importance has been 



the Haken-Strobl-Reineker model (Haken and Reineker (1972); Haken and Strobl (1973)) which in- 
corporates EVC by means of adding stochastic modulations to the bare electronic transitions energies 
and Coulomb couplings. This implies a certain model for the correlation functions of the energy gap 
coordinates, Eq. (12), (U m ,eg(t)U n , eg (0)), where the time-dependence is due to the dynamics of the 



nuclear degrees of freedom and the averaging is performed with respect to their equilibrium ensem- 
ble. More recently, these correlation functions became amenable to molecular dynamics simulations. 
Moreover, quantum-classical approaches have been developed to obtain the one-exciton Hamiltonian 



including the Coulomb coupling on-the-fly along classical molecular dynamics trajectories (Zhu et al. 
( |2007]|2008b|a| ) 



A different perspective is obtained by introducing collective harmonic normal mode vibrations 



( Renger et al. ( 2001 1 ; May and Kiihn ( 201 1 ) ) . For the purpose of illustration we consider the case where 
the dependence of the Coulomb coupling on the nuclear degrees of freedom can be neglected. The 
introduction of normal modes enables an analysis in terms of oscillators whose equilibrium positions are 
shifted upon electronic excitation, a problem which to some extent facilitates an analytical treatment 
(Huang-Rhys model, see Fig. [I]). In case of a continuous distribution of oscillator frequencies, EVC 
is described by means of a spectral density. Such a model applies in particular to the low-frequency 
part of the oscillator spectrum. Without discriminating between localized (at some monomer) and 
delocalized modes, the EVC Hamiltonian is given by ( |May and Kuhn] pOTTj )) 



Hi , y V 



E E h m + E E [<wg, 



m £ 



+ (l-S mn )(g^ e + g^} e )]q e \mn)(mri\ 



(23) 



Here {<?{} denotes the set of vibrational normal modes and gn/ a is the coupling constant for the £th 



mode and state \tp ne )- It follows from a linear expansion of the gap function, Eq. (12), in terms of q^. 

The absorption spectra of the dyes, which are used for making molecular aggregates are typi- 
cally characterized by a pronounced vibrational progression giving rise to a peak in the 1000-1500 
cm -1 range to the blue of the 0-0 transition (see Fig. Eh. The underlying coupling to higher fre- 
quency intramonomer modes often needs to be accounted for on a different footing. Here a possible 
strategy is the separation of such intramolecular modes, Q m ,(,, from the total set of aggregate normal 
modes which results in two EVC Hamiltonians plus some residual coupling between the (artificially) 
separated modes. The intramolecular modes can then be treated, for instance, non-perturbativcly 
by defining respective vibronic exciton states. While this is straightforward for a dimer or trimcr 
( |Matro and Cina| fl!995| ); |Kuhn et al] ( |1996a[ ); |Seibt et al.\ ( |2006| |2007| ) ; |Polyutov et al.\ ( |2011[ )), the 
treatment of larger aggregates requires further approximations concerning possible excitations. For 



instance, building on earlier work by Philpott (19711, Spano and coworkers employed a two-particle 



S 



approximation to a complete multi-particle basis to develop a Frenkel exciton polaron theory ( Spano 
(2010)). Here a local vibronic excitation is "dressed" by vibrationally excited ground state monomers 
An extension of the two-particle approximation has been proposed by Engel and coworkers using the 



multi-configuration time-dependent Hartree (MGTDH) ansatz for the correlated wave function (Seibt 
et al. (2009)). However, due to numerical limitations this approach seems to be restricted to smaller 
aggregates. Still a different strategy has been followed by Briggs et al. who employed a Green's func- 
tion approach in the so-called "coherent exciton scattering" (CES) approximation, which amounts to 



neglecting vibrational excitation in the electronic ground state (Roden et al. (2009b I ) 



2. 5. Dissipative EET Dynamics and Spectroscopy 

The EET dynamics in molecular aggregates is an example for a system with decoherence and 



dissipation (May and Kiihn (2011 )). The coherent dynamics of the (vibronic) excitons is disturbed by 



the interaction with some environment (heat bath) . Depending on the ratio between the Coulomb and 
system-bath couplings, the limits of (partly) coherent and incoherent transfer may be realized. The 
incoherent transfer limit under the assumption of strong intramonomer EVC is commonly known as 



Forster theory (Forster (1948)). A unifying framework for the theoretical description is provided by 
density matrix theory (Renger et al. (2001 ); May and Kiihn (2011 )). Here a reduced density operator, 
p{t), is defined by performing the trace of the total statistical operator with respect to the bath degrees 
of freedom. There are various methods for incorporating the system-bath interaction. Among the early 



approaches is the stochastic Haken-Strobl-Reineker model (Haken and Reineker (19721; Haken and 
Strobl|(p73|) which assumes a Gaussian-Markovian behavior of the bath-induced fluctuations. While 



this is a high- temperature theory, a generalization to finite temperatures has been given by Capek 
(Capek (1985)). More recently the non-Markovian case had been considered assuming a dichotomic 
behavior of the noise (Chen and Silbey (2011[)). The assumption of a stochastic influence of the 



environment is also present in the stochastic Schrodinger equation approach (Roden et al. (2009a)). 
While in most applications model assumptions concerning the fluctuations of transition energies and 



Coulomb couplings are introduced, they can be also obtained from quantum-classical simulations ( Zhu 
eTZ1 ( [2M7||2008a|b[ )) 



In the weak system-bath coupling limit multilevel Redfield theory formulated in terms of exciton 
eigenstates can be used, which is based on second-order perturbation theory assuming a Markovian 
dynamics (Renger et al. (2001 ); May and Kiihn (2011 )). Here, the bath enters via its spectral density 
which can be obtained, e.g., from fitting experimental spectra (Kiihn et al. (2002)). The approxima- 



tions of Redfield theory are lifted within the path integral (Thorwart et al. (2009)) and the hierarchy 



equation (Ishizaki and Tanimura (2005)) approaches, what comes at the expense of a restriction to 



small aggregates. 

To connect the quantum dynamics with spectroscopy, the macroscopic polarization needs to be cal- 
culated, which for a homogenous sample with aggregate density n agg is given by P(£) = rt agg tr(d ( o(t)), 
where the trace is with respect to the relevant system degrees of freedom and d is their dipole moment 
operator. In the weak field case one can expand the reduced density operator p(t) in powers of the 
external field to obtain the optical response in terms of linear and nonlinear response functions, i.e. 



(multi-)time correlation functions of the dipole moment operator (Mukamel (1995)). The response 
can also be expressed in terms of Green's functions, e.g. ( Chernyak et al. ( 1995 ); Mukamel and Abra- 



mavicms 



(2004); Roden et al. (2008)) or sum-over-states expressions (Knoester and Spano (1996)) 



A related approach is the so-called nonlinear exciton equation of motion method (Axt and Mukamel 



(1998); Mukamel and Abramavicius (2004)). Here a set of Heisenberg equations of motion for exci- 
ton creation and annihilation operators is derived and numerically solved in the time-domain. The 
advantage is that one obtains a clear picture concerning the elementary excitation contributing to a 
particular nonlinear signal. Further, being formulated in the local monomer basis, it scales favorably 
with aggregate size. The drawback is in the approximate description of dissipation and coherence 
and in the necessity to find appropriate truncation schemes for the many body hierarchy of equations 

While originally being 



of motion (Axt and Mukamel (1998); Mukamel and Abramavicius (2004)) 



formulated for electronic two-level monomers only, higher excited states (e.g. three-level model) can 
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be incorporated within the frame of the deformed boson theory (Kiihn et al. (1996b)). Increasing 
the field strength, the perturbative response function approach becomes cumbersome. In order to 
keep track of the field order and signal propagation direction, a carrier wave expansion of the density 
operator can be performed what leads to a hierarchy of equations for some auxiliary density operators 



(Renger et al. (2001)). 



2.6. Exciton-Exciton Annihilation 

It has been known for decades that annihilation of excitons in molecular aggregates is an important 
loss mechanism in EET. In luminescence experiments it leads to a reduction of the quantum yield 



and lifetime, thus restricting applications in optical devices. Early work by Suna ( 1970[ ) on molecular 
crystals focussed on the derivation of the bimolecular annihilation rate 7 in the equation for the 



exciton density n(r, t) given by (cf. related work by Kenkre (1980)) 



dn(r, t) 
dt 



-in 2 (r,t). 



(24) 



Subsequent work on linear aggregates by Malyshev and coworkers emphasized the role of higher excited 



electronic monomer states and golden rule type annihilation rates have been derived (Malyshev et al. 



(1999 2000 1). This problem has also been cast into a nonlinear exciton equation of motion approach 
based on a local site representation (Renger et al. (2000)). A consequent multi-exciton density matrix 
formulation of EEA has been given by May and coworkers ( |Renger et al. ( 2001 ); Bruggcmann and May 
(2003)). Invoking some approximations the two-exciton (ct^) to one-exciton (ai) interband relaxation 
rate can be cast into the simple form (Briiggemann and May (2003)) 



.(EEA) _\p 



n C l 2 £-( IC ) 



(25) 



(10 

Here, k m is the monomeric internal conversion rate for site to. Although analytic expressions for 



this rate can be given in terms of matrix elements of the non-adiabatic coupling operator, Eq. (13), 
any first principles calculation is hardly possible due to the fact that very highly excited electronic 
states will be involved. An important point of Eq. ( 25 ) is that the annihilation rate is proportional 

This does not come 
Recently, EEA has 



to the overlap of the local exciton wave function coefficients in the aggregate, 
unexpectedly in view of the local character of the internal conversion process 
been studied from the perspective of higher order perturbation theory of EET (May (20091). 



3. Applications 



1937) and Jelly (1936 1937) on 



3.1. Aggregation of Dye Molecules 

The first observations of aggregates were made by Scheibe ( 1936 
pseudoisocyanine dyes. Since then aggregation has been found to occur in a large number of molecular 
systems. Several excellent review articles provide overviews about the relevant compounds and their 
specific aggregation behavior (Mobius (1995); Bohn (1993); Wiirthner et al. (2011)). Cynanine dyes 



are probably the most popular group of molecules to study aggregates. Various cationic carbocya- 
nine (Scheblykin et al. (2001); van Burgel et al. ( |1994 ); Nemeth et al. (2009)), oxa- and thiacyanine 
( |Kuhn and Mobius| < |1971| )), ps eudoisocyanine (PIC) (|Sundstr5m et al.\ ( |1988[ ); |Fidder et al.\ ( |1991[ )), 
and neutral merocyanine dyes ( Hamanaka et aL| (2002 1) have been investigated. Chlorophyll dyes and 
structurally related compounds have been intensively studied due to their close relationship with nat- 



ural light harvesting systems (You et al. (2005); Balaban et al. (2005); Misawa and Kobayashi (1999J). 
During the last decade aggregates made from perylene bisimides (PBI) have attracted increasing in- 
terest, since these dyes exhibit high quantum yields and good photostability and are versatile building 



blocks for functional supramolecular structures (Kaiser et al. (2007); Li et al. (2008); Marciniak et al. 

poTTI)). 
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Figure 2: PBI-1 (7V,A r '-Di(A''-(2-aminoethyl)-benzamide)-l,6,7,12-tetra(4-tert-butylphenoxy)- 3,4:9, 10-perylene- 
biscarboximidc) monomer absorption spectrum and calculations for a reduced system: (a) leading HOMO-LUMO 
transition contributing to the So — > Si absorption, (b) absorption and fluorescence spectrum and calculated results 

narrow lines), (c) vibr ational mode that leads to the vibrational progression for the electronic So — > Si absorption 

Ambrosek et al.\ poTi) ) 
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Figure 3: PBI-1 monomer excited state absorption (ESA) spectrum, Si — > S n , which has been obtained by subtracting 
the ground state bleach (dashed grey) and stimulated emission (dotted grey) bands from the transient absorption (solid 
grey) spectrum measured at a delay time of 100 ps in dichloromethane (for more details, see |Ambrosek et al.\ fZ()U\). 



In Fig. [2] we show the analysis of the absorption spectrum of the perylene bisimide dye molecule 
A r ,A^'-Di(iV-(2-aminoethyl)-benzamide)-l,6,7,12-tetra(4-tert-butylphenoxy)-3,4:9,10-perylenebiscarbo- 
ximide (PBI-1) ( Ambrosek et aL| (2011 )), which includes the So — > Si as well as the So — > S2 transitions. 
The transitions can be characterized as being of tttt* type with the 7r electron density being located 
at the central PBI core with some contribution on the bay substituents. An important feature of 
the absorption spectrum is the shoulder on the blue side of the So — > Si transition which has been 
found to be due to a Franck-Condon progression. A normal mode vibration which contributes to this 
progression is shown in panel (c) of Fig. [2] Such signatures of electron- vibrational coupling are typical 
for many chromophores used for building molecular aggregates. It has strong implications also for the 
aggregate spectra and dynamics as will be shown below. 

In the context of nonlinear spectroscopy but also for the understanding of EEA the characterization 
of higher electronic states, S„, is important and in particular in the range of about twice the So — > S± 
transition energy. This spectral range is accessible by employing transient pump-probe spectroscopy 
(cf. Section 3.6 1 from which ESA profiles can be extracted. However, the actual procedure requires 
detailed knowledge of the ground state bleach and stimulated emission spectra for subtraction from 
the total signal. If this cannot be obtained spurious features will appear in the excited state absorption 
which facilitates a semi-quantitative discussion of the spectrum only. As an example we show the ESA 
of PBI-1 in Fig. [3j where at most two broad bands can be identified that extend from 500 and 600 
nm and from 650 to 750 nm. The theoretical assignment of these spectral features is hampered by the 
lack of accurate methods for highly excited states of these large organic molecules. For instance, using 
the DFT/MRCI approach, but for a reduced model of the PBI-1 monomer (axial groups replaced 



by hydrogens), transitions in this range are found at 861 nm, 665 nm, and 594 nm (Ambrosek et al. 

( poITT )). 

The attraction between the monomers of dye aggregates results from van-der-Waals interaction 
between the extended 7r-systems of the chromophores and from intermolecular hydrogen bonds. In 
particular the later are used to control the structure of the aggregates which is on the other side 
strongly influenced by steric effects. To obtain a detailed notion of the aggregate structure is a 
nontrivial task, but to some extent can be guided by Computational Chemistry. Exemplarily we show 
the result of a structural optimization of a PBI-1 hexamer in Fig. [4] 
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Figure 4: Structure of a PBI hexamer obtained from geometry optimization using the tight-binding DFT method. 
Clearly visible are the hydrogen bonds, which give some stabilization of the structure in addition to the van der Waals 
interaction of the it systems of the chromophores ||Ambrosek and Kiihn| ||2011|l). 



Many cyanine dyes seem to form a brick stone structure, at least when monolayers are prepared 
(Duschl et al. (1988); Mobius (1995)). For PIC aggregates investigations over several decades indicate 
that double stranded mesoaggregates are formed over which the excitons can be completely delocalized, 
at least at low temperatures (Daltrozzo et al. (1974)). The mesoaggregates can build up fibers which 



arrange in complex bundles at high concentrations (von Berlepsch et al. (20001; von Berlepsch and 



Bottcher (2002)). However, the exact structure depends on temperature, solvent properties, counter 



ions and probably also on the preparation procedure. 

3.2. Steady State Absorption and Fluorescence of Aggregates 

Aggregation is in most cases demonstrated by changes in the stationary absorption spectrum in 
comparison to the monomer in highly diluted solutions. In the case of formation of J-aggregates a 
sharp band, red shifted with respect to the first electronic absorption band of the monomer is observed 
which rises with increasing dye concentration. In J-aggregates the slip angle between the monomeric 
transition dipole moments and the intermolecular connection axis is below the magic angle of 54.7° 
(cf. 



Eq. (|5[)) and the lowest exciton state carries most of the oscillator strength (Bohn ( 1993 ); 
Burgel et al. (1994)). Figure [5] shows the changes in the case of a PBI aggregate (Li et al. (2008); 
Marciniak et al. (2011)). The lowest band gains additional absorption strength, experiences a slight 



red shift, and becomes narrower. Here the derealization length is only about two monomers (see 
below). For large coherence domains very sharp absorption bands can result ( |Mobius| ( [T995| ). 

H-aggregates exhibit broad absorption bands which are blue shifted with respect to the monomer 
absorption since here the slip angle is larger than the magic angle and the upper excitonic states are 



responsible for the absorption ( |Bohn (1993); Wiirthner et al. (2011)). 

In Fig. [6] we have plotted results on the energy level structure and oscillator strengths distribution 
for a vibronic dimer as a function of the Coulomb coupling strength ( |Polyutov et al. ( 2011 )). The model 
includes the possibility of excitation of one vibrational quantum in the electronic excited state of each 



13 



24000 22000 20000 18000 



16000 



50000 - 



o 
E 



wavenurnber (crn 1 ) f} 
/ ; 


i 


0.16 mM / \ 




0.0016 mM / ; 




/,-< 


k I 


t 


\\ 




\\ 


, g -|i, */ 






\\ 









450 500 550 600 650 
wavelength (nm) 



Figure 5: Spectral changes due to the formation of J-aggregates of a core-substituted perylene bisimide dye | |Li et aL\ 
( |2008} ; [Marciniak et aT] ( |2011| l). 



monomer. For small coupling strengths (J/w v ib) J/H-aggregates have most oscillator strength in the 
transition to the lowest/highest state corresponding to the 0-0 transition. For J-aggregates this holds 
true also for increasing J/uj vi ^, whereas for H-aggregates the exciton- vibrational state which carries 
oscillator strength is strongly mixed with respect to electronic and vibrational contributions. This 
simple model already hints at the more complex absorption line shape for H-aggregates as compared 
to J-aggregates (Eisfeld and Briggs (2006); Spano (2010)). 

In several cases a splitting of the J-band in two or more components is observed and attributed 



to the coexistence of two or more locally slightly different arrangements (Cooper (1970); Fidder 
( 1991| )). If the structure of aggregates consists of unit cells containing two molecules on 



et al. 



non- equivalent sites Davydov splitting occurs. In the case of 3,3'-bis(sulfopropyl)-5,5'-dichloro-9- 
cthylthiacarbocyanine (THIATS) this results in the simultaneous appearance of J- and H-bands in 
the absorption. Excitation of the higher lying H-band leads to J-type fluorescence demonstrating that 
intcrband relaxation from the H- to the J-states takes place (Scheblykin et al. (2000)). 

The linewidth of the J-band is strongly reduced compared to the monomer absorption (see above) . 
This results from mo tional narrowing caused by the delocalization of the excitons over many monomers 
(Knapp et al. (1984)). The width is roughly TV -1 / 2 times the monomer width where N is the number 
of monomers within the coherence domain. Applying this relation to PIC aggregates a coherence 
domain of about 60 monomers was estimated (Knapp et al. (1984)). The domain size is often also 
called coherence length since the electronically excited states of chromophores within such a domain 
are coherently coupled. Numerical simulations of the absorption band taking energetic disorder into 
account indicate an exciton delocalization over about 100 molecules ( Fidder et aL| ( 1991 1 ) . 

Another elegant way to characterize the number of coherently coupled chromophores is saturation 
spectroscopy with short laser pulses. From the saturation fluence 7 sa t = hv/a one can determine 
the absorption cross section a of the absorbing entity and can compare it to the cross section of 
the monomer. The short pulse duration ensures that the ground state is not repopulated during the 
excitation. In this way the coherence domain was estimated to be about 100 monomers for 1,1'- 
diethyl-2,2'-cyanine chloride (PIC-C1) in water (Kopainsky and Kaiser (1982)) and about 340 for 
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Figure 6: Upper panel: Energy level structure of the one-exciton states, a\ = 1, . . . , 4, for a vibronic dimer as a function 
of the Coulomb coupling. The arrows indicate possible intraband transitions as discussed in Fig. |9| L ower panels: 
Oscillator strengths for one-exciton transitions in H- and J-aggregate configuration i jPolyutov et al\ ( |201 1 [ ) ) . 



l,l'-diethyl-2,2'-cyanine bromide (PIC-Br) in a mixture of water and ethyleneglycol (Kobayashi and 



Sasaki (1994|). 



The fluorescence originates from the lowest exciton state since intraband relaxation is very efficient 
(see below). Accordingly, the fluorescence spectrum is a mirror image of the J- absorption band 
(Mobius ( 1995)). In an ideal J-aggregate the Stokes shift should be zero since the lowest exciton state 
carries most of the oscillator strength and motional narrowing suppresses the impact of optical active 
vibrations. Energetic disorder leads in fact to a distribution of optical accessible excitonic states. 
However, the width of this distribution is typically narrow and in many cases extremely small Stokes 
shifts are observed like e.g. in PIC aggregates (Fidder et al. (1991)). In the case of weakly fluorescent 



dye molecules like PIC-C1 the formation of J-aggregates can drastically increase the fluorescence yield 
(Katheder (1940)). This results from the increase of the transition dipole due to coupling of many 
chromophores and probably also from the reduced flexibility of the monomer within the aggregate 
which can suppress internal conversion pathways. The excitons are mobile and can change between 
different derealization domains. This results in an overall relaxation towards the lowest lying exciton 
states and can cause a significant Stokes shift of the fluorescence (Scheblykin et al. (2001)). The 
mobility of the excitons can be nicely demonstrated by incorporating small amounts of molecules in 
the aggregate which function as energy acceptors for the excitons. In this case extremely efficient 
quenching of the aggregate fluorescence occurs. It was found in PIC ag gregates that one quencher 
molecule can deactivate 10 3 to 10 6 chromophores (Scheibe et al. (1939)) exceeding by far the size 
of a coherence domain. This shows that an exciton can visit several domains within its lifetime. In 
bleaching experiments applying optical near-field microscopy an upper limit for the exciton diffusion 
length of 50 nm was found at room temperature for aggregates made from the dye 1, l'-diethy 1-2,2 '- 
cyanine iodide (PIC-I) (Higgins and Barbara (1995)). Incorporation of electron accepting or donating 



molecules in an aggregate can also efficiently quench the aggregate fluorescence (Penner and Mobius 



(1982)). If both kinds of molecules are incorporated electron transfer from the donor to the acceptor 
can occur via the aggregate. In this way it was demonstrated that charges can migrate up to 4 nm in 
aggregated monolayers of PIC dyes. 



15 



The fluorescence of H-aggregates is usually weak or totally lacking ( Herz ( 1977 1 ) . This results from 
the fact that the optically accessible exciton states are in the upper region of the exciton band and 
their population is quickly transferred to low-lying exciton states which possess no transition dipole to 



the ground state (Fink et al. (2002)). Interestingly, excitons in nanowires made from H-aggregates of 
perylene tetracaboxylic diimide seem to exhibit a reasonable mobility. In these wires defects imprinted 
by a scanning microscope tip lead to effective quenching which points to an exciton diffusion length 



of at least 250 nm (Chaudhuri et al. (2011)) 



3.3. Exciton Lifetime 

The lifetime of excitons in aggregates depends on various factors. As mentioned above the forma- 
tion of J-aggregates can suppress nonradiative decay channels of the monomer. However, the radiative 
lifetime of J-aggregates is strongly reduced compared to the monomer since the coherent coupling of 
many chromophores results in a large transition dipole moment and the radiative rate is proportional 
to the number of chromophores in a coherence domain (Bohn ( 1993[) ; van Burgel et al. (1994)). This 
effect is often called superradiance (Knoester and Agranovich (2003)). The increase of the radiative 
rate due to aggregation was used to estimate the coherence length of zinc chlorin aggregates which 



are models for light harvesting systems to 10-15 chromophores (Prokhorenko et al. (2002)). 

The radiative rates of H-aggregates are very low since the transition of the lowest exciton state 
to the ground state is dipole forbidden (see above). The exciton lifetimes are usually limited by 
nonradiative processes. H-aggregates of perylene tetracaboxylic diimide show a lifetime increase from 



1.7 ns to 46 ns upon cooling from room temperature to 25 K (Chaudhuri et al. (2011 )). This indicates 



that the exciton decay is caused by a thermally activated nonradiative process. 

3.4- Temperature Dependence 

The temperature influences the excitonic properties mostly in two ways. Firstly, with increasing 
temperature geometric and energetic fluctuations and phonon scattering decrease the size of coherence 
domains and reduce the number of effectively coupled monomers. Secondly, if the temperature is very 
low hops between different domains are suppressed. Then, in case of significant static disorder, the 



excitons can not reach the lowest excitonic states but are trapped at local minima (Scheblykin et al. 



(|200_Lj)). 

In the case of Langmuir-Blodgett films made from l-methyl-l'-octa-decyl-2,2'-cyanine perchlorate 



the fluorescence lifetime r changes from 8.2 ps at room temperature to 5.5 ps at 143 K ( Dorn and Muller 
(1987)). The data were interpreted applying the equation 1/r = A: ra d(l + f3T) ( |Mobius and Kuhn 
(1979)). Here, fc ra d is the radiative rate of the aggregate and T the temperature. (3 is a characteristic 
parameter of the particular dye, which parameterizes the nonradiative decay rate fc nr = fc rac j/3T and 
describes its reduction with decreasing temperature. The data evaluation indicates that the size of 
the coherence domains increases by a factor of about 3 upon cooling from room temperature to 143 K 
( |Dorn and Muller| ( |1987| |). 

For J-aggregates of THIATS a minimum of the Stokes shift at intermediate temperatures of about 



110 K and a complex temperature dependence of the fluorescence bandwidth were found ( jScheblykin 
et al. (2001)). The relative Stokes shift which is the Stokes shift divided by the bandwidth of the 
fluorescence shows an interesting behavior. It is constant at very low temperatures up to 30 K, decays 
then exponentially with temperature, and adopts at about 80 K a pretty temperature independent 
value again. This was interpreted as an indication that at low temperatures the excitons are trapped 
at local energetic minima (|Scheblykin et al. (2001)). At intermediate temperatures the hopping 



probability and the mobility of the excitons increases more and more with temperature till global 
minima can be reached. A similar behavior was also observed in core-substituted PBI aggregates 
(Kaiser et al. ( 2009[ )). There the critical temperatures were a factor of 2.8 higher than for THIATS 
pointing to a higher degree of disorder. The nonmonotonic temperature behavior of the Stokes shift 
due to the interplay between local bandstructure determined by the static disorder and relaxation 



dynamics has also been investigated theoretically by Bednarz et al. (2003). 



16 




Figure 7: Dependence of the two-photon absorption (TPA) spectrum of a heterodimer (transition energy difference is 
0.05 eV, J = —0.05 eV) composed of electronic three-level systems on the local anharmonicity (A n ) and the ratio of 
transition dipole moments (k„). (a) K n = 0.5, (b) A n = 0, (c) A n = 0.05 eV (after Kiihn et al. (1996b)). 



3. 5. Two-Photon Spectroscopy of the Two-Exciton Band 

The two-exciton band is characterized by an interplay between local and nonlocal double exci- 
tations. Their mixing depends on two parameters, i.e. the local anharmonicity of the electronic 
transitions, A„ = huj n e f — HLu n ^ eg , and the ratio between transition dipole moments between So — > 
Si and Si — > S n transitions, K n (assuming that they have the same orientation). In Fig. [7] we 
show calculated two-photon absorption spectra of the two-exciton band for a heterodimer com- 
posed of three-level monomers. Increasing the anharmonicity parameter (panel (a)) the eigenstates 
forming the two-exciton band become less mixed and one can identify the delocalized double ex- 
citation around u)\ + cj 2 = 3 eV and the two local double excitation at higher energies. Increas- 
ing the ratio K n leads to an increased mixing between the eigenstates causing level splitting since 
Jmn oc d m j e d n ,eg = Kmd m ,egd n ,eg ■ Finally we notice that the nonlinear response vanishes if the 
system is harmonic, i.e. when A n = and n n — v2 (panel (b)). 



3.6. Transient Spectra 

Ultrafast pump-probe spectroscopy is very helpful to characterize excitonic properties of aggre- 
gates. Figure [8] shows a transient absorption spectrum obtained 10 ps after optical excitation of 
core-substituted PBI aggregates (Marciniak et al. (2011)). The strong negative feature is due to the 



bleach of the ground state absorption and stimulated emission from the electronically excited state 
(see broken lines in Fig. [8| . Its sharp blue wing results from a partially hidden photoinduced absorp- 
tion. It can be extracted by subtracting the bleach and the stimulated emission from the transient 
spectrum. The obtained ESA (grey line) shows a dominant band which is in shape and intensity 
similar to the low energy band of the ground state absorption but slightly blue shifted. This is the 
signature of a one- to two-exciton transition. The corresponding ESA band in aggregates of PIC-Br 



was directly observed at 1.5 K due to the very narrow bandwidth ( jFidder et al. (1993)). The weaker 
bands at shorter wavelengths in Figj8] result from the vibronic structure of this transition while the 
absorption band at longer wavelengths is due to an ESA that is also observed in transient absorption 
measurements of the monomer (cf. Fig. [3]). 

The difference, A.E, between the transition energies of the ground state absorption and the tran- 
sient absorption band caused by the one- to two- exciton transition can be approximated for linear 



aggregates by (van Burgel et al. (1994)) 



AE: 



3ir 2 J 

(N + iy 



(26) 
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Figure 8: Transient absorption spectrum (black line) 10 ps after optical excitation of PBI aggregates and its decompo- 
sition in ground state bleach, stimulated emission, and excited state absorption. 



Corresponding simulations of the transient absorption spectrum for aggregates made from 5, 5', 6,6'- 
tctrachloro-l,l'-diethyl-3,3'-di(4-sulfobutyl)benzimidazolocarbocyanine (TDBC) reveal at room tem- 
perature a delocalization length of about N — 15 monomers (van Burgel et al. ( 1994 1 ) . At a tem- 
perature of 1.7 K a delocalization length of 50-100 was estimated from transient absorption data for 
PIC-Br aggregates (Durrant et al. ( 1994[ )). For a merocyanine dye which forms two dimensional aggre- 
gates the coherence domain at 4.2 K was characterized by a combined analysis of the stationary and 
the transient absorption spectrum (Hamanaka et al. (2002)). The results indicate that the domain 



extends over 26 molecules in one direction and over 40 in the other direction. 

Finally, we point out that effects of coherent exciton motion and dephasing can also be detected 



as distinct features in ultrafast pump-probe experiments (Heijs et al. (2007)) 



3. 7. Photon Echo Experiments 

Photon echo experiments are suitable to investigate line broadening mechanisms acting on the 
J-band. In the temperature range from 1.5 K to 100 K the photon echo decay measured at PIC-Br 



aggregates accelerates from 14 ps to about 4 ps ( De Boer et al.\ (1987)). This was attributed to 



dephasing caused by the coupling to a low-frequency mode. However, even at 1.5 K the photon echo 



decay is faster than the fluorescence lifetime (Fidder et al. ( 1991 )). It was argued that several excitonic 
states contribute and interference effects reduce the photon echo signal in time. The decay of the two- 
pulse photon echo signal of TDBC aggregates at room temperature was analyzed simultaneously 
with the stationary J-absorption band using a fluctuation correlation function with two components 
described by the magnitudes (A) and time scales (A -1 ) (van Burgel et al. (1994)) 



(u n , eg (t)u n , eg (o)) = A; 



low 



-A sl , 



vt i a 2 — Af aBt t 
T ^fast e 



(27) 



The fast component, Af ast , was found to be in the range of 5 fs with an amplitude of 10 THz and 
the slower one, A s i ow , in the range of 10 ps with an amplitude of 54 THz. Overall, the dephasing 
is slower than for dye monomers in solution indicating that the delocalization of the excitons over 
several chromophores averages out the fluctuations. 



3.8. Intraband Relaxation 

Each excitonic band is a dense manifold of excitonic states characterized in case of linear aggregates 
by their quasi- momentum k. Higher lying fc-states can be populated by photo-excitation. Relaxation 
to lower fc-states is quite efficient due to the small energy differences. Two-dimensional spectroscopy is 
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Figure 9: Population relaxation within the vibronic one-exciton manifold provided that the highest state is initially 
prepared (cf. arrows in Fig. [61. (a) J/u v -ib = 2 (curves at 1 ps correspond to states ai = 4, 3, 2 with increasing value), 
(b) J/aj v ib = 0.3 (curves at 1 ps correspond to states a\ = 2, 3, 4, 1 with increasin g value). The upper panel shows the 
effect of transient population trapping (for more details, see |Polyutov et al.\\201l)). 



particular suited to study intraband relaxation since it provides directly the excitation dependence of 
the dynamics. Using this technique it was found that in the case of l,l'-diethyl-3,3'-bis(sulfopropyl)- 
5,5',6,6'-tetrachlorobcnzimidacarbocyanine aggregates the rate for intraband relaxation accelerates 



from zero at the red wing of the absorption band to l/(40fs) in the blue wing (Stiopkin et al. (2006)). 

The intraband relaxation might be influenced by the vibronic substructure of the exciton band. As 
an example the population relaxation in the one-exciton band according to Redfield theory is shown 
in Fig. [9] for the vibronic dimer model of Fig. [6j Most notable is the so-called transient population 
trapping observed in the case of strong Coulomb coupling (panel (a)). With increasing ratio J/w v ib 
the level separation increases (cf. Fig. [6]) such that the transition frequency to reach the lowest state 
moves out of the frequency range where the spectral density of the coupled bath modes differs from 
zero. Thus on the time scale given, no equilibrium population is established within the one-exciton 
vibronic manifold (Polyutov et al. (2011)). Finally, we point out that the effect of vibronic coupling 
can also be viewed in real space, i.e. in the context of exciton transport. Here it has been reported 



that vibronic coupling may have a strong influence on the exciton coherence size (Rodcn et al. (2009b); 



Womick and Moran (2011)) 



3.9. Interband Relaxation 

The two-exciton band consists of states in which two excitons are on the coherence domain and 
are coupled to monomeric double excitations. Up to now the population of a two-exciton state in 
aggregates has not yet been unambiguously observed. This is usually attributed to the short lifetime 
of such a state. A blue shifted component of the transient absorption features was found during a short 



time interval after the excitation in the case of merocyanine J-aggregates (Hamanaka et al. (2002 1 ) 



and of PIC-Br aggregates (Minoshima et al. (1994)). It was attributed to the annihilation induced 
decay of multiple excitons occupying one coherence domain. However, at a very low temperature of 
1.5 K the width of the one- to two-exciton transition in PIC-Br aggregates is only 5 cm -1 (Fidder et 



al. (1993)). This indicates that two-exciton states with a lifetime of at least 1 ps should exist under 
these conditions. 

The microscopic origin of the interband transition is the breakdown of the Born-Oppenheimer 
approximation. According to Section[2.6|it is the local character of the exciton states which determines 
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Figure 10: Population dynamics of the two-exciton states of the heterodimer model of Figs. [7^, (A n = 0.2 eV, panels 
(a)) andJTJ; («„ = 2, panels (b)) which follows from a hierarchy equation of motion approach including an internal 
conversion rate of 1/500 fs. Notice that due to intraband relaxation in the case of panel (b) there is some intermediate 
population trapping in state 0.2 = 1 which has a pronounced nonlocal character (|Yan and Kuhn| (|2011[)) . 



the transition rate, Eq. (25). To demonstrate this in Fig. 10 we show the population dynamics of 
the heterodimer model of Fig. [7] as obtained using a hierarchy equation of motion approach ( Yan 
and Kuhn (2011)). For simplicity it is assumed that the highest state of the two-exciton band is 
populated. Panels (a) correspond to the case of Fig. [7^, (A„ = 0.2 cV), i.e. local and nonlocal double 
excitations are weakly mixed. Since the initial state has a pronounced local character the decay via 
interband transitions is rather rapid. In contrast in panel (b) the case of Fig. [Tj; (n n — 2) is used, 
where nonlocal and local states are strongly mixed. Having a more pronounced nonlocal character 
the initial state decays much slower. 

Intraband relaxation is also very important if several excitonic J-bands exist within the same aggre- 
gate. In the case of porphyrin J-aggregates a decay time of about 1 ps was observed for the relaxation 
from the exciton band formed by the molecular S2-states down to the band resulting from the molec- 
ular Si-states ( |Misawa and Kobayashi (1999)). J-aggregates of Tetrakis(4-sulfonatophenyl)porphyrin 
(TPPS4) exhibit also two J-bands, a strong one around 490 nm and a weaker one at 707 nm. The 
exciton coherence length in the upper band is about 10 chromophorcs. After optical excitation of 
this band ultrafast decay to the lower band occurs and the coherence length reduces to one or two 
molecules (Gulbinas et al. ( 2007| )). Subsequently the exciton relaxes to a not or only weekly radiating 
state within 200 ps in agreement with the observation that the fluorescence yield decreases with aggre- 
gation. This kind of self-trapping is not very common in J-aggregates. In TPPS4 the transition dipole 
to the lowest excited electronic state is small and thereby also the excitonic coupling what might favor 
localization of the exciton. In complex aggregates with several overlapping J-bands two-dimensional 
spectroscopy is a suitable tool to disentangle various contributions to the relaxation dynamics as well 
as interband coherences (Milota et al. (2009)). For an aggregate with four overlapping bands a hi- 



erarchy of relaxation steps in the range of 100 fs to 1 ps was observed in this way (Nemcth et al. 



(2009)). 



3.10. Exciton- Exciton Annihilation 

Laser spectroscopic experiments found that the exciton lifetime in J-aggregates decreases with exci- 
tation density and time resolved measurements revealed that the associated kinetics is non-exponential 
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Figure 11: Transient absorption signal (grey dots) of excitons on PBI J-aggregates for different excitation energies. The 
observed kinetics is modeled by assuming EEA and one- (black solid line) and three-dimensional diffusion (broken line) 
of the excitons (see also |Marciniak et al\ pOllJ ) . 



(Sundstrom et al. (19881). These phenomena occur already at moderate excitation densities and re- 
sult from EEA. Due to their mobility excitons can approach each other and interact. In such an 
encounter one exciton is deactivated with a high probability while the other remains. The mechanism 
was already introduced in the Introduction and in Section 2.6 In the simplest case the time depen- 
dent exciton density n(t) obeys the rate equation, Eq. (24 1, supplemented by a term describing the 
intrinsic decay with rate k\ . 

For a signal S which is proportional to n plotting dS/dt/S versus the signal itself should give a 
linear rise. For PIC chloride a deviation from this behavior was observed at low excitation densities 
indicating that annihilation occurs only above a certain threshold density (Sundstrom et al. (1988)). 
From this density one can estimate the size of the domain over which an exciton migrates and in 
which at least two excitons have to be generated to induce annihilation. In this way it was shown 
that in PIC chloride excitons can sample up to 10 4 monomers (Sundstrom et al. ( 1988)). For THIATS 
exciton migration over 6 x 10 4 monomers at room temperature and even over 6 x 10 6 at 77 K was 
observed ([Scheblykin eta!] (|1998[)). 



A more sophisticated analysis was carried out for EEA on PBI J-aggregates (Marciniak et al. 
(2011)). Figure 11 shows the time dependent absorption signal due to optical excitation of PBI 
aggregates in methylcyclohexane for different excitation energies. The negative signal is proportional 
to the number of excitons and is obtained by integrating the spectral region from 590 nm to 625 nm, 
where bleach and stimulated emission are dominating the transient absorption spectrum (cf. Fig. 
[8]). The exciton signal exhibits a non-exponential decay which accelerates with increasing excitation 
density indicating that EEA takes place. The time-dependent exciton concentration is modeled taking 
annihilation into account and assuming that the excitons can diffuse within the aggregates. In this 
case 7(i) of the above rate equation becomes time dependent and adopts the form 7id(0 — y/2D/nt 
for one-dimensional and 73d(0 = 4:TrRD(l + R/ ^/2nDt) for three-dimensional diffusion. D is the 
diffusion constant and R the critical distance between two excitons when annihilation takes place. For 
the three-dimensional case, insufficient agreement with the data is achieved while the one-dimensional 
model reproduces the measurements almost perfectly (see Fig. 111. This indicates that the excitons 
move along the PBI aggr egates as on one-dimen sional strings. The analysis reveals also a diffusion 
constant of 1.29 nm 2 /ps (Marciniak et al. (2011 1). It demonstrates that annihilation dynamics can 
be used to investigate local transport properties. 
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4. Concluding Remarks 

Excitation energy transfer in molecular aggregates has a long history, but only recent years have 
witnessed an unprecedented step forward in our mechanistic understanding of the underlying pho- 
tophysical properties. Ultrafast nonlinear spectroscopy provides insight into the real time dynamics 
of exciton motion, unravelling its coherences and their decoherence due to the interaction between 
excitonic and nuclear degrees of freedom. For example, in the past the concept of an exciton coherence 
domain has been inferred from indirect indicators such as band shapes or emission lifetimes. Nowadays 
two-dimensional spectroscopy provides a much more direct look and is capable to disentangle time 
scales and correlations of fluctuations. This development comes along with a much increased prospect 
for realistic modeling. Still the Frenkel exciton Hamiltonian and models like the dimer or the linear 
chain are indispensable. However, the development in Computational Chemistry allows for the first 
time to obtain first principles based parametrizations for concrete systems. Last but not least, organic 
synthesis has taken up the challenge to design and synthesize dye aggregates with desired electronic 
properties. The convergence of these different areas holds great promises for the future development, 
which includes, for instance, the integration of molecular aggregates into organic solar cell devices. 
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